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ABSTRACT 

We investigate the nature of the diffuse Galactic emission in the Wilkinson 
Microwave Anisotropy Probe (WMAP) temperature anisotropy data. Substan- 
tial dust-correlated emission is observed at all WMAP frequencies, far exceed- 
ing the expected thermal dust emission in the lowest frequency channels (23, 
33, 41 GHz). The WMAP team (Bennett et al. ) interpret this emission as 
dust-correlated synchrotron radiation, attributing the correlation to the natural 
association of relativistic electrons produced by SNae with massive star forma- 
tion in dusty clouds, and deriving an upper limit of 5% on the contribution of 
Draine & Lazarian spinning dust at K-band (23 GHz). We pursue an alternative 
interpretation that much, perhaps most, of the dust-correlated emission at these 
frequencies is indeed spinning dust, and explore the spectral dependence on en- 
vironment by considering a few specific objects as well as the full sky average. 
Models similar to Draine & Lazarian spinning dust provide a good fit to the 
full-sky data. The full-sky fit also requires a significant component with free-free 
spectrum uncorrelated with Ha, possibly hot (~ lO^K) gas within 30° of the 
Galactic center. 

Subject headings: cosmic microwave background — diffuse radiation — dust, 
extinction — ISM: clouds — radiation mechanisms: thermal — radio continuum: 
ISM 



1. INTRODUCTION 

The first year of data from the Wilkinson Microwave Anisotropy Probe (WMAP) has 
sparked a revolution in the study of cosmology and the early universe (Bennett et al. 2003, 
Hinshaw et al. 2003, Spergel et al. 2003). In addition, it has provided full sky maps of the 
interstellar medium at 23 — 94 GHz at a sensitivity of ~ 200/iK per 7' pixel at resolutions 
much higher than previous full-sky observations. The quality of the data has rekindled 



^Hubble Fellow, Henry Norris Russell Fellow 



- 2 - 



the debate about the origin of "anomalous" dust-correlated microwave emission (Kogut et 
al. 1996, dc Oliveira-Costa 1997,1998,1999,2002,2003; Finkbeiner et al. 2002, 2003; Banday 
et al. 2003; Lagache 2003) dubbed "Foreground X" by de OUveira-Costa et al. (2002). 

The WMAP team relied on a Maximum Entropy Method (MEM) analysis to determine 
that spinning dust^ with the Draine & Lazarian (1998) cold neutral medium (CNM) spectrum 
accounts for less than 5% of the emission in any WMAP waveband. We show that for 
different models of spinning dust, these constraints do not apply. Indeed, in most of the 
Galaxy, spinning dust emission (or some other mechanism such as magnetic dipole emission, 
Draine & Lazarian 1999) explain the data just as well. 

We emphasize that the WMAP team's primary objective was to prevent Galactic fore- 
ground emission from contaminating their cosmological results based on the CMB anisotropy. 
They have done a superb job of modeling the foregrounds, regardless of whether the phys- 
ical interpretation of the subtracted temperature anisotropy signal is correct. The WMAP 
polarization data also provide important cosmological constraints as well as additional infor- 
mation about foregrounds. For example, if the Bennett et al. hard synchrotron component 
can be shown to be ~ 30% polarized, then spinning dust would become an unlikely explana- 
tion (see Lazarian & Draine 2000). On the other hand, if it is unpolarized, the synchrotron 
would be ruled out. Something between the two extremes is likely, and this will help us 
determine the relative contributions of synchrotron and spinning dust. As the polarization 
data are not yet public, discussion of polarization will be deferred to a later paper. 



2. ANALYSIS 

A disadvantage of the MEM analysis used by the WMAP team is that it can only find 
the components it is searching for. Significant spectral variation is expected not only in the 
synchrotron spectrum, but also in any spinning dust that might be present. The WMAP 
strategy, using data at the WMAP frequencies and priors based on the Haslam et al. (1982) 
408 MHz map, the Finkbeiner (2003) Ha map, and the Finkbeiner et al. (1999; hereafter 
FDS99) dust prediction, produces a "synchrotron" template highly correlated with thermal 
dust emission (see Fig. 4 in Bennett et al. 2003). The explanation for this is that the 
most energetic electrons are found nearest dusty star-forming regions, so there is a natural 



^The idea of grain rotational emission was first discussed by Erickson (1957). More recently, after the 
discovery of the population of ultra small grains, Ferrara & Dettmar (1994) noted that the rotational emission 
from such grains may be observable, but their treatment incorrectly assumed Brownian thermal rotation of 
grains. 
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correlation between dust and synchrotron emission. This is certainly plausible, but it is not 
the only possible explanation for the strong correlation with dust. 

Instead, we start by asserting that there are three components whose spectral charac- 
teristics are well enough understand to use external data sets, subtract these off in fits with 
very little freedom, and then approach the remaining emission with an open mind. These 
three components are free-free, soft synchrotron, and thermal dust. Here we differentiate 
between "soft" synchrotron (traced by the Haslam map at 408 MHz) and "hard" synchrotron 
as described by Bennett et al. (2003). The existence of the soft synchrotron is not in dispute. 



Although warm ionized gas is well traced by Ha emission in three recent surveys (VTSS, 
Dennison et al. 1998; SHASSA, Gaustad et al. 2001; WHAM, Reynolds et al. 2003) as 
presented by Finkbcincr (2003) and corrected for dust using the Schlegel, Finkbeiner & Davis 
(SFD; 1998) dust map, one must worry about several systematics. Bennett et al. (2003) 
point out several sources of uncertainty, paraphrased as follows: (a) Ha light scatters off 
dust, while microwave free-free scattering is negligible; (b) Ha calibration is uncertain at 
the 10% level; (c) for the 1/4 of the sky with no WHAM coverage, there is significant 
uncertainty (~ IR) separating geocoronal Ha from Galactic ISM emission; (d) The assumed 
Gase B Balmer atomic rates may not be universally valid; (e) The ratio of Ha to free-free 
microwave emission depends on Tg, which may vary from object to object (e.g. Heiles et 
al. 2000); (f) There is uncertainty in the ionization state of He; (g) ionized gas with bulk 
velocity outside of the WHAM velocity coverage is not included in the Finkbeiner map; and 
(h) the extinction correction assumes that the ionized gas and dust arc co-extensive along 
the line of sight. Of these sources of error, (e) and (h) likely dominate, while (b), (c), and 
(f) are well characterized, and (a) and (d) are probably small. 

In spite of these difficulties, we can still make use of an Ha template. The first step is 
to limit the Ha extinction predicted by the SFD (1998) dust map to be A < 2 mag. For dust 
with optical depth r = A/(1.086 mag) uniformly mixed with Ha emitting gas the observed 
intensity is 



In the limit of large r, the intensity is simply reduced by a factor of r. For A = 2 mag the 
uniform mixing case results in a reduction factor of 2.2, between the extremes of dust in 
front (factor of 6.8) and dust behind gas (factor of unity). Though uniform mixing appears 
to be a reasonable approximation in general, we shall bear in mind that factor of ~ 2 errors 
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in the free-free correction are likely in some parts of the sky. 

A systematic variation in Tg is also expected, with lower temperatures toward the Galac- 
tic center, where higher metallicity allows for higher cooling efficiency. Lacking an electron 
temperature map, we assume Tg = 5500 K everywhere, and use the free-free "haze" (§3.3) 
to absorb the errors in this. Using the conversion formula in Valls-Gabaud (1998), an EM 
of 1 cm-6 pc corresponds to 0.391 r-i oi72o-o.o29/T4 Rayleigh, (7^ = 7;/10^ K) which equals 
0.636 R for T = 5500 K. 

Once an emission measure (EM) is derived for a given pixel on the sky, the spectrum 
is constrained by physics (Spitzer 1978) so only one free parameter (Tg) is needed for the 
whole sky. This is valid between the plasma frequency cutoff and the optically thick limit 
at low frequencies. According to catalogs (Kurtz, Churchwell, & Wood 1994) of optically 
thick ultra-compact H ii regions the filling factor of such regions is so small that at WMAP 
frequencies a negligible fraction of the free-free is expected to be optically thick, and certainly 
within our rather conservative mask there can be no significant ultra-compact Hii regions. 

2.2. Soft Synchrotron 

Synchrotron emission from relativistic electrons with energy distribution N{E) ~ E''^ 
produces a synchrotron flux density spectral index a = —(7 — l)/2. In temperature units 
the index is f3 = — (7 -|- 3)/2. 

The electron energy distribution is expected to vary spatially in the Galaxy, with a 
harder spectrum near the sources of energetic electrons, e.g. supernova remnants, and a 
softer spectrum elsewhere. Bennett et al. (2003) argue that the natural correlation of SNR 
with dusty star-forming clouds should cause a strong correlation between hard synchrotron 
emission (at ~ 30GHz) and dust, with the correlation lessening at lower frequencies, so that 
at 408 MHz, the Haslam (1982) "soft" synchrotron map shows very little correlation to dust, 
other than the usual concentration of emission in the Galactic plane. This is certainly a 
cogent argument that synchrotron emission should become more dust-correlated at higher 
frequencies, but provides no evidence that such radiation dominates other mechanisms, such 
as spinning dust. 

There is little doubt, however, that soft synchrotron (traced by the 408 MHz survey) 
is present, and important at least at K band (23 GHz) so we include it in the removal of 
"known components". A brightness temperature index oi fi — —3.05 removes the most 
obvious K-band synchrotron features at high latitude, so we assume this index is constant 
throughout the high-latitude sky. 
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2.3. Thermal Dust 

Another component whose morphology and (approximate) spectral characteristics are 
understood is thermal (vibrational) dust emission from grains large enough to be in equilib- 
rium with the interstellar radiation field. Emission from this dust peaks at ~ 140/im and 
deviates strongly from a thermal blackbody spectrum. A Rayleigh- Jeans emissivity function 
of has often been assumed in the hterature (e.g. Draine & Lee 1984; Schlegel, Finkbeiner 
& Davis 1998) but when dust temperature variation is accounted for, the COBE FIRAS 
data (Fixsen et al. 1998) are better fit by a steeper power law emissivity (/? = 2.6) near 
the peak and = 1.7 at lower frequencies, with a break at about 500 GHz (Finkbeiner et 
al. 1999; FDS99). This fit ties the IRAS and DIRBE data to FIRAS via a fit with only 4 
global parameters describing the two emissivity laws, and the requirement that the emission 
is dominated by grains in equilibrium with the interstellar radiation field. Predictions at 6' 
resolution based on this fit are available on the web.^ 

Because the spectral slope used to extrapolate the FDS99 thermal dust spectrum through 
the WMAP bands is relevant to this work, we briefly review its justification. The FDS99 
two-component model yields a substantially better fit (reduced = 1-85 compared to 31 for 
a u'^ model) even when the spatial and spectral covariance of the FIRAS data are included 
(FIRAS Explan. Supp., Brodd et al. 1997). The FIRAS data cannot differentiate between 
a single component with non-power-law emissivity and two independent power law compo- 
nents, but the physical reality of the two separate components is not implausible. Amorphous 
silicates with a wide range of emissivity indices /? ~ 1.2 — 2.7 have been observed in the lab 
(Agladze et al. 1996), including amorphous MgO-2Si02 which has a very high microwave 
emissivity to optical absorption ratio, leading to rather different mean temperatures (9K and 
16K) for the two components. These lab emissivities were measured at ~ 300 GHz and 20K 
and may become steeper at lower frequencies. However, this interpretation of the spectral 
break is hardly unique; if the dominant emitter has such a break in its emissivity function at 
500 GHz, then a single component could explain the data just as well. Another explanation 
that has been advanced is very cold dust grains spatially mixed with the warm dust - (Reach 
et al. 1995), though a physical mechanism for keeping the grains so cold is not provided. 
Such a model would predict a steeper slope at lower frequencies also. 

Regardless of interpretation, the Finkbeiner et al. (1999) model has been very successful 
in the sub-mm - microwave, though small but interesting deviations from the model have been 
observed by BOOMERANG (150, 240, 410 GHz; Masi et al. 2001). At lower frequencies, 
however, there is a surprise. 
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Comparing FDS99 predictions to COBE DMR, Finkbeiner et al. found that COBE 90 
GHz was slightly higher, but at 53 and 31 GHz the emission per dust column is a factor of 
2.2 and 31 higher than expected. These results are similar to the earlier Kogut et al. (1996) 
results derived without an explicit dust temperature correction. Because of this it was 
expected that the FIRAS-bascd predictions would agree well with WMAP 94 GHz, but 
be significantly contaminated by some other dust-correlated emission mechanism at lower 
frequencies. This appears to be true. 

Even though we can expect the vibrational thermal dust brightness temperature to 
go as T ~ with j3 = 2 through the WMAP channels, there may be other emissivity 
enhancements such as magnetic dipole emission. These are treated as a separate component, 
so that temperature and other environmental dependence can be searched for also. 

The WMAP team (Bennett et al. 2003) has determined that the spectral slope from 94 
to 61 GHz is /3 = 2.2, significantly steeper than the j3 = 1.7 extrapolated from the Finkbeiner 
et al. (1999) fit. However we find in this work that /3 = 1.7 is consistent with the data, and 
lacking any good theoretical reason to depart from the FDS parameterization, we continue 
to use this index in the following. 

2.4. Spinning dust template 

For a spinning dust template we use an empirically determined power of dust color tem- 
perature (derived from the DIRBE 100/240/im ratio) times column density. This template 
refiects the observed Tdust dependence of the anomalous emission, but assumes constant spec- 
tral dependence across the sky. We use T|^^^ times the FDS prediction at 94 GHz, which is 
proportional to dust column times Tj^^^ because the Rayleigh-Jeans dust prediction already 
contains one power of temperature. This modulation by varies the dust template by a 
factor of ~ 2 peak-to-pcak for the SFD temperature range (16 — 20K). There is no obvious 
theoretical reason for this temperature dependence other than the vague notion that hotter 
dust tends to be in a more exciting radiation field, or perhaps has a size distribution with 
more small grains, either one of which might enhance spinning dust emission. We will fit 
this average spectrum to the (nearly) full sky and then examine a few regions of interest to 
see deviations from it. 
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2.5. CMB anisotropy 

The CMB dipole is subtracted by the WMAP team, and used to cahbrate the anisotropy 
pattern to < 1% uncertainty. The WMAP sensitivity is high enough that CMB anisotropy is 
clearly visible in all bands and (for this work!) represents a significant nuisance background 
signal. Because the cosmological signal is not correlated with the foreground morphology for 
any physical reason, the cross correlation between the foreground templates and the CMB 
(due only to chance alignment of features) is minimal. Nevertheless, we use the superposition 
of the 5 bands suggested by Bennett et al. (2003) to estimate the CMB anisotropy, and 
subtract it from each band. The superposition used is 

ATcMB = -lOQii: - .684X0 - .096Q + 1.921^ - .250W^ (2) 

where {K, Ka, Q, V, W) refer to the thermodynamic AT in the 5 WMAP bands respectively. 
This works quite well except where bright Galactic point sources near the ecliptic plane 
show the beam asymmetries.^ This superposition only imperfectly cancels the foregrounds, 
of course, allowing them to contaminate the CMB estimate and leak back into our fore- 
ground analysis. We find that any such contamination is negligible, as it would appear as a 
foreground with the spectrum of the CMB, and no such residual is found in our analysis. 

2.6. Mask 

We mask out the 208 sources listed in the WMAP source catalog, and also the LMC, 
SMC, M31, Orion-Barnard's Loop, and NGC 5090 (a total of 2.7% of the sky). Additionally, 
areas where the Ha extinction given by SFD is greater than 2 mag, i.e. A(Hq;) = 2.65E{B — 
V) > 2, are masked (8.9%). The mask is smoothed to 1° FWHM to match the smoothed 
data, excluding a total of 11.4% of the sky from the analysis 

2.7. The Fit 

The fit minimizes in 8-dimcnsional parameter space over unmasked pixels. First 
all data and templates are smoothed to 1° resolution using spherical harmonic convolution 



*It is often said that the inherently asymmetric WMAP beams are symmetrized by scanning each point 
on the sky from many different approach angles. This is true near the ecliptic poles, but at low ecliptic 
latitude the beam asymmetry is obvious in the K and Ka band maps. 
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on the N — 512 HEALPix^ sphere. Noise is derived from ctq values and A^obs, the number of 
observations in each pixel (Jarosik et al. 2003). Initial guesses arc determined by extensive 
experimentation with the data, then passed to the simple-minded but fairly robust amoeba 
routine from Numerical Recipes (Press et al. 1992) as implemented in IDL.® Although the fit 
is a simultaneous regression on all parameters, it is instructive to "peel away" the foregrounds 
one at a time and show each WMAP band at each step. 

Figure 1 is a 6 column by 7 row grid of images, each a Cartesian projection of half 
the sky (-90 < / < 90°; -90 < b < 90°). The first five columns are the WMAP bands in 
frequency order (K, Ka, Q, V, W) = (23, 33, 41, 61, 94 GHz). The sixth column is the 
template being subtracted from the current row to yield the images on the row below. The 
first column is stretched so that -1 mK is black and 1 mK is white (= ±15 Jy/sr). Columns 
2-4 are stretched to the same intensity limits (in Jy/sr) as col. 1, and col. 5 matches the 
stretch of column 4 in thermodynamic temperature units (to avoid saturation of the gray 
scale). Therefore, the CMB appears to become stronger from columns 1 through 4, but is 
the same brightness in col. 5 as col. 4. 

Row 1: the sky signal (after CMB dipole removal) for the 5 WMAP bands, with the 
plane and point sources masked. Col. 6 is the CMB map described above. 

Row 2: CMB-subtracted sky; note the Hii region around COph, above the plane and 
just left of center. This feature is prominent in the Ha based free-free template, shown at 
right. The free-free spectrum is constrained by physics, so only one parameter is fit for an 
overall normalization. 

Row 3: Sky minus CMB and Ha; note the "Loop I" synchrotron feature in the north, 
visible in the Haslam-based soft synchrotron template at right. A single power law is fit 
from the Haslam survey frequency (408 MHz) to WMAP. The spectrum is likely concave 
downward over such a large range in frequency, but WMAP 23 GHz dominates the fit, and 
is the only frequency for which this component is important anyway. 

Row 4: Sky minus CMB, Ha, and soft synchrotron. Dust correlated emission is evident, 
as traced by the Finkbeiner et al. (FDS; 1999) template at right. This template is subtracted 
with amplitudes determined by the FDS99 fit to FIRAS, not by cross-correlation with the 
WMAP data. Note that the subtraction at 94 GHz appears to be nearly perfect. 



'''sec http:/ /www.eso.org/scicncc/healpix. Smoothing is done with spherical harmonic convolution on the 
HEALPix sphere, using smoothing code implemented in IDL by D. Finkbeiner. 

^IDL is a product of Research Systems, Inc. http://rsinc.com 
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Row 5: Sky minus CMB, Ha, soft synchrotron, and FDS dust. There is still signifi- 
cant dust-correlated emission present at lower frequencies. The FDS dust emission template 
times T^„5^ (SFD dust temperature squared, based on DIRBE 100 and 240/im emission ratio) 
appears to be a good template for the remaining component, tentatively interpreted as spin- 
ning dust, and is shown at right. A cross- correlation coefficient is determined independently 
for each WMAP band. 

Row 6: Residual, after subtracting CMB, Ha, soft synchrotron, FDS dust, and the 
putative spinning dust, is nicknamed "the Haze." It looks nearly the same in every band 
and has the spectrum of free-free emission. A template for the Haze is produced from the 
average fit residuals in K and Ka bands, and subtracted from all five WMAP bands with a 
free-free spectrum. There are no free parameters in this operation, other than one value per 
sky pixel determined from the K and Ka maps. 

Row 7: The WMAP residuals after subtracting all above templates show very little 
remaining other than measurement noise. The one exception is residual emission from COph 
which appears to peak at 41 GHz and roll off to either side (relative to free-free). The point 
source mask used in the analysis is shown in the lower right hand panel. 

3. INTERPRETATION 

The three important results of this analysis are that (1) the quality of the fit is very 
high, (2) the dust-correlated component has a reasonable spectrum to be spinning dust, and 
(3) a "haze" unassociated with Ha appears to be present in the data. 

3.1. Quality of fit 

The quality of the fit is remarkable. All that is subtracted from each sky map is a 
CMB template; a map with free-free spectrum (the Ha based template plus the haze); 
soft synchrotron, which is hardly controversial at low frequencies and unimportant at high 
frequencies; FDS dust emission extrapolated from the FIRAS fit; and emission strictly cor- 
related with FDS dust times T^y^^^. This last component is the only one that can contain 
hard synchrotron emission or spinning dust emission - the others are all well constrained 
and physically understood. The surprisingly small residuals imply that the physical mecha- 
nism primarily responsible for this emission is closely tied to the dust grains themselves (e.g. 
spinning dust) rather than a mechanism relying on the spatial correlation of dust and some 
other emission source (e.g. hard synchrotron emission). 
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in doing this fit, we have not allowed the spectral dependence of any component to 
float from pixel to pixel, as the WMAP team MEM analysis did (Bennett et al. 2003). 
Though it may seem overly conservative to fix the spectral dependence of each component, 
this rigidity in the fit allowed the discovery of the free-free haze, and the quality of the 
resulting fit indicates that spectral variation of the components is not required to fit WMAP 
data off the plane. Because the templates were applied directly rather than as priors to a 
MEM analysis, evaluation of is straightforward. Values given in Table 1 are for 1° FWHM 
gaussian beams; this reduces measurement noise to approximately 15/iK per beam in each 
band. Because of the CMB and haze subtraction, the residual in each band is actually 
a superposition of the foreground templates and 5 WMAP bands appropriately weighted. 
The smoothed WMAP noise for each band is then combined in a weighted quadrature sum 
to produce the expected noise in the 5 residual maps. This process assumes that the 5 
WMAP bands have uncorrelated noise, and that the Ha, Haslam, and FDS templates have 
no measurement noise. The fact that per dof is of order unity at high latitude (Table 1) 
indicates that these assumptions are justified in Q, V, and W bands, and less so in K and 
Ka bands. At \h\ < 30° the x^/dof is formally poor (13.6 at K and 7.8 at Ka, and 2 at Q) yet 
represents removal of 98%, 97%, and 99% of the variance in the WMAP data respectively. It 
is significant that these values are with respect to the 1° FWHM smoothed data; the residual 
variance of 60//K per beam at K band is less than half of the median measurement noise per 
pixel, and relatively smaller in the other bands. 

The residuals in the Bennett et al. (2003) analysis are also impressively small. It is 
impossible to compare the results of this study to the Bennett et al. MEM analysis, which 
fit several numbers for each pixel on the sky (constrained by prior templates). Such a MEM 
analysis does not allow for an easy interpretation of x^, so a smaller fit residual would not 
imply the superiority of one foreground parameterization over the other. In fact, the same 
freedom that allows the MEM to minimize foreground contamination of the cosmological 
signal may actually inhibit the physical interpretation of the foreground signal by allowing 
an overly simple model to fit the data. 

3.2. Spinning Dust Spectrum 

Draine & Lazarian (1998, 1999) chose to phrase the comparison between theory and 
observation in terms of emissivity per H atom, using units of Jy cm^ sr~^ per H in their 
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work^. This ratio is appealing to theorists because it directly relates the dust model to 
chemical abundances in the ISM. However, measuring the hydrogen column density, includ- 
ing molecular, atomic, and ionized forms, is technically difficult. Finkbeiner et al. (2002) 
measured cross-correlation slopes of emission against the SFD98 dust map, and converted to 
iV(H) using 1 mag E{B — l^) = 8 x 10^^ ^(H)- Although this conversion is appropriate for 
diffuse high latitude material, the gas to dust ratio is known to vary in general, and phrasing 
measurements in terms ol E{B — V) or A'"(H) when they have nothing to do with either has 
led to confusion. 

In this work, we stay as close as possible to the empirical data and state emissivity per 
dust in terms of optical depth at 100/xm, or rioo/i- This optical depth is defined with respect 
to a blackbody at the reference dust temperature T^e/ = 18.175 K used in SFD98. This 
temperature is the median high-latitude dust color temperature determined from the DIRBE 
100//m / 240//m intensity ratio assuming a emissivity law. For the SFD calibrations, tioo^x 
of unity corresponds to 266 mag E[B — V) or 2.13 x 10^^ -^(H), so the reader may convert 
Draine & Lazarian models and previously published measurements to ju/Tioof^ in Jy sr~^ by 
multiplying by 2.13 x 10^^. 

The cross-correlation coefficients for the FDS99 xT^^^^ spinning dust template are given 
in Table 2. These may be added to the FDS99 thermal dust prediction, where the sum is 
taken at dust temperature Tref = 18.175 K. In Figure 2 the dust-correlated emission per Tiqo^ 
is shown for 5 Draine and Lazarian (1998) models of electric dipolc emission from spinning 
dust, along with the cross-correlation points from WMAP and other projects. Note that the 
thermal dust emission and "spinning dust" emission have been combined for WMAP, just 
as they have been for other experiments (intentionally or not). The WMAP data points 
fall about a factor of two lower than the WNM model. These WMAP points represent the 
average high-latitude sky, so there is no reason to expect them to agree with any particular 
Draine & Lazarian model, but rather a superposition. 

Two spinning dust candidates previously studied with the Green Bank 140 foot tele- 
scope (Finkbeiner et al. 2002) are also shown, Lynds Dark Nebula 1622 (Figure 3, see also 
Hildebrand & Kirby 2003) and LPH 201.663+1.643 from the Lockman, Pisano, & Howard 
(1996) catalog of diffuse Hii regions (Figure 4). A superposition of the WNM and CNM 
models plus a small amount of free-free yields a reasonable fit for LDN1622. Note that the 
thermal dust level is not being fit, merely extrapolated from the far IR according to the 
FDS99 FIRAS fit. Spinning dust models have a few free parameters to tune, and some 



''Plots in the early Draine & Lazarian articles are incorrectly labeled "Jy sr ^ per H" , an error repeated 
by Finkbeiner et al. (2002). 
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models may come closer to fitting the LPH cloud, though this will require further theoretical 
work. Lacking a physical model for this object, a curve (the Draine & Lazarian WNM model 
shifted down a factor of 1.8 in frequency and up a factor of 3.7 in amplitude) is shown in 
Figure 4 for comparison. Neither of these objects would be well fit by a superposition of 
synchrotron, free-free, and thermal dust without including spinning dust. 

3.3. The Haze 

A rather substantial amount of emission (both positive and negative) was attributed 
above to the "free-free haze," without any interpretation. As the spectrum of this component 
most resembles free-free (Figure 5), the straightforward interpretation is that it constitutes 
the error in the Ha based free-free template. Because the spatial template is orthogonal 
to the others by construction, it contains positive and negative parts. Some error is to be 
expected both from variation in warm gas temperature about the mean Tg^s — 5500K used, 
and from the fact that some gas is hot enough (> lO^K) that the Ha recombination fine 
is strongly suppressed. This need for hot gas applies especially to the region south of the 
Galactic center, where there is significant emission seen in the haze but none of the other 
templates have a similarly extended feature of significant brightness. 

Further study of this diffuse haze, perhaps using O vi absorption measurements toward 
distant stars below the Galactic Center, may be able to verify this hypothesis about its 
origin. We emphasize that the tentative identification of the haze as free-free is based only 
on its spectral shape, and if WMAP polarization data show a significant polarization, that 
would suggest very hard synchrotron (~ t/~°'^) emission as a source. Further consideration 
of the haze is deferred to another paper. 

A HEALPix map of the K-band haze template used in Figure 1, defined as {Tk + 
2.186 Tko)/'^-! is available on the web^. 

4. CONCLUSIONS 

We have performed an analysis of the WMAP temperature anisotropy data, using free- 
free, soft synchrotron, FDS thermal dust, and spinning dust as templates. We find: 

1. The FDS thermal (non-magnetic) dust prediction is good on average over the high- 



'http://skymaps.info 
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latitude sky at 94 GHz, and even in dark cloud LDN1622 and H ii region LPH 201.663+1.643. 

2. Free-free is fit by T = 5500K ionized gas (rather cool, but see e.g. Heiles et al. 2000), 
plus a haze toward to the Galactic center with free-free spectrum. This is the first detection 
of emission from this component. Further study is needed to test whether this component 
is indeed free-free or perhaps a very hard synchrotron component. 

3. The residual dust-correlated emission is well fit by T^^^^x dust column density over 
most of the sky. Such a model, along with free-free, soft synchrotron, and FDS thermal dust, 
accounts for ~ 98% of the variance in the WMAP data. 

4. There is no confiict between the general idea of spinning dust emission and the 
WMAP data. In spite of the Bennett et al. (2003) statement that the Draine & Lazarian 
CNM model accounts for less than 5% of the ISM emission (which is true), spinning dust 
remains a reasonable hypothesis both morphologically and spectrally. 

5. The Green Bank Galactic Plane Survey GPA data at 8 and 14 GHz (Langston et 
al. 2000) rule out hard synchrotron as a majority component in the diffuse ISM in the 
WMAP bands (see companion paper, Finkbeiner, Langston, & Minter 2003). This leaves 
spinning dust as a viable alternative. 

The WMAP data alone cannot establish the relative levels of spinning dust, magnetic 
dust, and hard synchrotron emission, since any of these theories can be adjusted to explain 
the data. However, objects such as LDN1622 (Figure 3) are clearly not dominated by 
synchrotron in the WMAP bands, and are well fit by a superposition of spinning dust models. 
Therefore, the dominance of spinning dust at 23 GHz is a good working hypothesis. Further 
data (including surveys at lower frequencies, better WMAP signal to noise, and WMAP 
polarization data coming soon) will help lift the remaining degeneracies and determine the 
identity of the dust-correlated microwave emission from"Foreground(s) X." 

I am indebted to Bruce Draine, Ed Jenkins, and David Schlegel for encouragement and 
helpful advice. Roger Hildebrand called the WMAP spectrum of LDN1622 to my attention, 
and provided other insights. Lyman Page, Chuck Bennett, and Gary Hinshaw provided 
helpful comments on the proper use of the WMAP data in this work. Comments from 
an anonymous referee substantially improved the text. Some of the results in this paper 
were derived using HEALPix (Gorski, Hivon, and Wandelt 1999). This research made use 
of the NASA Astrophysics Data System (ADS) and the IDL Astronomy User's Library at 
Goddard^ DPF is a Hubble Fellow supported by HST-HF-00129.01-A. 



^ http : / /idlastr o . gsf c . nasa. gov / 
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Fig. 1.- 



The WMAP foreground grid; see detailed discussion in §2.7. 
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Fig. 2.— Model dust emissivity per rioo^ for DC, MC, CNM, WNM, and WIM conditions 
(as in Draine & Lazarian 1998b, Figure 9) with {solid lines) and without {dashed lines) 
contribution from vibrational dust at mean temperature. For comparison with the Draine & 
Lazarian models, unit tioo^ corresponds to A^(H) = 2.13 x 10^"^ cm~^. Gray line is emission 
from free-free for < UeUp > / < uh >= 0.01 cm~^ averaged along the line of sight. WMAP 
data from Table 2 (including thermal dust) are overplotted {filled circles) Also shown are 
measurements from the COBE/DMIi {solid diamonds) from Finkbeiner et al. (1999), sim- 
ilar to Kogut et al. (1996); Saskatoon {open diamonds) (de Oliveira-Costa et al. 1997); 
the Cottingham & Boughn 19.2 GHz survey {open square) (de Oliveira-Costa et al. 1998), 
OVRO data {solid squares) (Leitch et al. 1997); Tenerife data {triangles) {de Oliveira-Costa 
et al. 1999); GB 140 foot (crosses) (Finkbeiner et al. 2002). The OVRO points have been 
lowered a factor of 3 relative to Draine & Lazarian (1998b, Figure 9), because the unusual 
dust temperature near the NCP caused an underestimate of rioo^ along those lines of sight. 
Given the large range of model curves, all measurements are consistent with some superpo- 
sition of spinning dust, vibrational dust, and free-free emission. 
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V (GHz) 

Fig. 3. — LDN1622. Green Bank (5, 8, 10 GHz) measurements in a 6' beam (Finkbeiner 
et al. 2002) and WMAP (23-94 GHz) data points smoothed to a 60' beam (symbols), each 
with statistical la error bars only. Thermal dust estimate from Finkbeiner et al. (1999) 
(dashed line), free-free estimate (dotted line), spinning dust (thin line), and sum of free- free, 
thermal dust, and spinning dust (thick line). The spinning dust model is a superposition of 
two Draine & Lazarian models: 0.6x WNM plus 0.11 x CNM model. 
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V (GHz) 

Fig. 4.— LPH 201.663+1.643. Green Bank (5, 8, 10 GHz) and WMAP (23-94 GHz) data 
points (symbols), thermal dust estimate from Finkbeiner et al. (1999) {dashed line), free- 
free estimate {dotted line), spinning dust {thin line), and sum of free-free, thermal dust, and 
spinning dust {thick line). None of the published spinning dust models come close to fitting 
this object. For reference, the CNM model is shown, shifted in frequency and amplitude. This 
curve may not represent a physical spinning model for any reasonable choice of parameters. 
Green Bank points are 6' and WMAP is smoothed to 60', so deviation of the fit from WMAP 
points may be due to beam mis-match between Green Bank and WMAP, but the outher at 
10 GHz is significant. 
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Fig. 5. — The Haze is determined in 4 WMAP bands by subtracting CMB, soft synchrotron 
(Haslam template), free- free (Ha template) and spinning dust. Using the K-band haze as a 
template, it is then subtracted from Ka, Q, and V bands assuming various power laws. A 
free-free spectrum fits most of the sky well, apart from the ( Oph cloud (/, 6) = (5°, 25°). 
^3.3. 



-22- 



Table 1. Chi squared results, single band 



map 


b cut 


K 


Ka 


Q 


V 


W 


Mean 




deg 


23 GHz 


33 GHz 


41 GHz 


61 GHz 


94 GHz 






(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


W 


all 


391.8645 


139.7677 


93.6773 


43.6718 


31.5367 


140.1036 


W 


< 30 


583.3032 


206.6596 


133.5638 


54.6387 


38.4774 


203.3285 


W-C 


all 


122.9073 


17.7629 


8.7606 


7.0739 


2.2803 


31.7570 


W-C 


< 30 


185.5578 


27.3719 


13.7034 


10.6546 


2.8084 


48.0192 


W-C-T 


all 


14.4250 


2.8392 


1.8906 


1.7544 


1.1566 


4.4132 


W-C-T 


< 30 


27.1182 


4.6733 


2.7873 


2.5053 


1.2265 


7.6621 


W-C-T-H 


all 


9.4000 


5.4792 


1.6167 


1.4913 


1.1282 


3.8231 


W-C-T-H 


< 30 


13.6420 


7.8132 


2.0239 


1.8904 


1.1632 


5.3065 


W-C-T-H 


> 45 


2.0515 


1.5306 


1.0962 


1.1211 


1.0938 


1.3787 



Note. — Col. (1): W=WMAP data, C=CMB map defined in Eq. 2, T=Template 
fits, H=Haze defined in §3.3, Col. (2): Galactic latitude cut, in addition to mask 
defined in §2.6, Col. (3-7): per dof for the 5 WMAP bands in 1° beams. Col. 
(8): mean for the 5 bands. In W band the is dominated by imperfect CMB 
removal. In other bands imperfect foreground removal is dominant, but the model 
reduces \b\ < 30° excess variance (above measurement noise) in the data by factors of 
46, 30, 130, and 60 in K-V bands respectively. 
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Table 2. Correlation slopes 



Band ' Freq. Spinning Dust Spin + Thermal 
GHz MJy sr-Vfoi^ MJy sr-VfoL 



(1) 


(2) 


(3) 


(4) 


K 


22.5 


4.33391 


4.37490 


Ka 


32.7 


3.18627 


3.34376 


Q 


40.6 


2.45289 


2.79613 


V 


60.7 


1.39604 


2.85618 


W 


93.1 


0.88224 


7.69205 



Note. — Col. (1): WMAP band name, Col. (2): 
effective central frequency for free-free spectrum, Col. 
(3): cross-correlation slope of spinning dust compo- 
nent per 100//m optical depth. Col. (4): sum of spin- 
ning dust and thermal (vibrational) dust predicted by 
Finkbeiner et al. (1999) per tiooh- These data are 
plotted in Figure 2. 



